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SUMMARY 


Kubo’s formulation of irreversible quantum statis- 
tics is applied to electron mobility. Xn the weak coup- 
ling limit the approach to equilibrium is shown to be 
characterized,, in the lowest order, by a relaxation time. 
With the use of Coopersmith wave functions for hard-core 
interactions a characteristic relaxation time is derived 
which to "all orders" is independent of the external 
(small) driving force. For the specific application to 
mobility the general theory yields the known results of 
semi-classical kinetic theory at high temperature and 
low density. For the region of lower temperature and 
higher density the theory gives qualitative agreement 
with the anomalous drop in mobility found by Levine and 
Sanders for electrons in low temperature helium. How- 
ever, lack of quantitative agreement indicates a need 
for further work on the problem. 


ABSTRACT 


Kubo's formulation of irreversible quantum statis- 
tics is applied to electron mobility. In the weak coup- 
ling limit the approach to equilibrium is shown to be 
characterized, in the lowest order, by a relaxation time. 
With the use of Coopersmith wave functions for hard-core 
interactions a characteristic relaxation time is derived 
which to "all orders" is independent of the external 
(small) driving force. For the specific application to 
mobility the general theory yields the known results of 
semi-classical kinetic theory at high temperature and 
low density. For the region of lower temperature and 
higher density the theory gives qualitative agreement 
with the anomalous drop in mobility found by Levine and 

Sanders for electrons in low temperature helium. How- 
ever, lack of quantitative agreement indicates a need 
for further work on the problem. 
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INTRODUCTION 


1 * 

It has long been realized that the quantum mechan- 
ical analogue to the Boltzmunn equation contained two 
assumptions of limited and even questionable validity . 1 
These are the repeated random phase approximation (RPA) 
and the relaxation time assumption (RTA) Recently sev- 
eral authors have worked at developing a general trans- 
port theory that would avoid these assumptions. Van 
Hove chose to start with the derivation of the master 
(Pauli) equation and showed that RPA could, for various 
physically reasonable potentials, be replaced by a much 

9 

weaker condition. A number of attempts have also been 
made starting from the Liouville equation governing the 
density matrix (> . Greenwood derived an equation for 
conductivity in a metal which avoided the use of RPA.** 
This formula was evaluated by Edwards for a metal with 
randomly located impurities to yield the usual solution 
to the Boltzmann equation.^ Lax also developed a formal 
theory which avoids use of RPA. For the case of a weak 
scattering perturbation (wealc-coupling limit) and with- 
out RTA he obtained the usual transport reiult and the 

5 

Nyquist theorem. Kohn and Luttinger have also shown 
the validity of the quantum mechanical Boltzmann Eq. 
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in the weak coupling limit without resorting to RPA 
or RTA. They carried their work to third order in the 
coupling strength parameter, and for a periodic lattices 
with impurities restricted to lattice sites found devia- 
tions from standard theory which appear to play a role 
in the Hall effect.** 

Kubo also treated irreversible processes starting 

7 

from the Liouville equation. He developed a formal 
theory using neither RPA nor RTA. Kubo»s theory bears 
much the same relationship to Irreversible Statistical 
Mechanics as does the partition function to equilibrium 
theory. It is an exact, formal theory with (to my know- 
ledge) no questionable assumptions and quite often is 
rather difficult to evaluate. This theory has been 
applied by Nakano to conductivity in metals who showed 

that it gives the Gruneissen formula if one introduces 
8 

RTA. Chester and Thellung have taken Kubo's formal 

Q 

theory and applied it to elastic scattering in a metal. 
For the weak coupling limit they obtained results ident- 
ical with those of the quantum Boltzmann theory with 
neither RPA nor RTA. 

In the- present work we calculate electron mobility 
based on the formal theory developed by Kubo. Neither 
RPA nor RTA are used. The weak coupling limit is con- 
sidered briefly in lowest order. Our primary concern, 
however, is the calculation of electron mobility to "all 
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orders" in the presence of randomly located harc-core 
scatterers. We shall see that in both instances the 
approach to a steady state is actually characterized 
by a single relaxation time which we derive, not assume* 

Our choice of a model based on random hard- sphere 
scatterers corresponds most closely to the situation in 
helium vapor. For low energy electrons, helium has a 
positive scattering length, no Ramsauer effect and the 
electron-atom interaction is adequately characterized by 
a spherically symmetric hard-core repulsion. 

Electron mobility in helium gas has been experimen- 
tally investigated at various times since the early 
1920s at temperatures ranging from 77° K to 300 K. 

The results of these experiments have been effectively 
explained by classical kinetic theory. Recently, there 
has been interest in electron mobility in low temperature 
helium gas. This was motivated by measurements in liquid 
helium where it was found that an electron formed a sta- 
ble dense complex with very low mobility. This led to a 
prediction that such a complex would also be stable in 
helium gas at sufficiently low temperature and high den« 
sity. Measurements to test this hypothesis were made by 
Levine and Sanders at temperatures near 4° K and densi- 
ties of the order of lCr atoms/cm . Under these con- 
ditions the mobility was found to be almost four orders 
of magnitude lower than that obtained from classical 


kinetic theory of free electrons. A theory based on 
formation of a correlated (bubble) state was developed 
for this high density- low temperature phenomenon, but it 
failed to account for the transition region at lower den*" 
sity. 

A quantum-mechanical consideration of the properties 

of a free electron experiencing hard-core repulsion in an 

ideal gas was made by Coopersmith in a calculation of the 

13 

equilibrium free energy. By using Kubo’s quantum mech- 
anical formulation for mobility, we extend this work to 

the study of a non-equilibrium property. The first order 

14 

results of this effort have been previously reported 
and are reviewed in this work in the course of develop- 
ment of the full theory. 

All work will be done in units with 'f) = m - 1 . 
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II. MOBILITY - GENERAL THEORY 
The mobility of an electron in the direction 
of an external electric field E is given b y^ 

fJ. = Ct (Mj. Af) ( 2 . 1 ) 

where tr is a trace and U/y is the electron velocity 
along E. A f is the difference between the total 
density matrix and the unperturbed (E=0) density 

matrix y 5 , 

Af - f* T - f 3 ( 2 - 2 ) 

The Hamiltonian H,j, for an electron in the pres- 
ence of N scattering centers can be correspondingly 
divided into 


H r = H + AH 

(2.3) 

where 



(2.4) 

AH--- )eE 

(2.5) 

and 



N.V . ^ 
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Vj,H V(% )-V(r r Pe) (2.6) 

As mentioned in the first section, Kubo started 
from the Liouville equation 

[Ht.fr] ( 2 . 7 ) 

and for a weak external perturbation obtained the line- 
arized adiabatic solution 


Af* -iZ 


-i 


*7 n . -‘fH -pH it" 

J n . eft q [AHjC ]e 


( 2 . 8 ) 


where t^ is the length of time the system has been 
acted on by A H and the equilibrium density matrix 
is given by 


Z exp 

(2. 9. a) 

z = th exp(-/3H) 

(2.9.b) 


As expressed in equation (2.8) Aj* is an average 
(integral) over the length of time the system is 
acted bn by the perturbing force. The averaging is 
done to the Heisenberg representation of the commutator 
of the perturbation and the equilibrium density matrix. 
A derivation of equation (2.8) is given in Appendix A. 
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Combining equations (2.1) to (2.9) gives 

<fte />,<? Je JJ)? ( 2 . 10 ) 

This equation will be the starting point for the pres- 
ent work. 

Kubo actually extended the formal development 

of the theory by interchanging the trace and the time 

integral. The resulting integrand (including the 

trace) he identified as the response function. This 

response function can then be treated as an average, 

with the trace defining the averaging process. It is 

this response function that Kubo and others have used 

2 

as the starting point for their calculations. While 
this procedure has calculational advantages for their 
work (e.g, cyclic permutation of operators), it does 
not lend itself to the present work. One consequence 
of taking the trace first is the elimination of infor- 
mation regarding the specific functional form of the 
time dependence. (Specifically, in the course of cal- 
culating we would like to consider the limits of 

validity of RTA. Since we expect that the relaxation 
time T when it does exist is momentum dependent, 
is best formulated in a way that displays this depend- 
ence explicitly. We have sufficient freedom to achieve 
this since the trace makes the choice of representation 


ss ie Z ' tk £ 


arbitrary. In light of the above comments the obvious 
procedure is to work in the momentum representation 
and retain the trace for the last step in the calcula- 
tion. Some simplification accrues as a result of the 
above choice, since is diagonal in the momentum 

representation. This follows from the definition of 
and equations (2.3) to (2.5). 

-^5 = ‘-[r'A, >1 r ^ (2. 11) 

In the momentum representation the eigenvalue equation 
is 

M> rili> (2. 12. a) 

and the normalization is 

<A /,?'>= Sl-A-1') , <h*> = (2. 12. b) 

where p is the momentum operator, k its corresponding 
eigenvalue, S(-A-Jl') the Direc delta function and S2 
the volume of the system. 

The above formulation would yield fA valid only 
for a particular specific arrangement of the scattering 
centers, with explicit dependence on the N position 
parameters y. } j,s /, 2, /V . The mobility we actu- 

ally want is the ensemble (configuration) average over 
all possible arrangements. For a system of N random 
(uncorrelated.) scatterers 


3 


9 


ju = JT cfpf 


(2.13) 


where the bar indicates an snsemble average. It will 
be assumed that the ensemble average commutes with the 
trace and the integration over t. 

In the representation in which p is diagonal, 
the matrix element arising from equation (3.10) 

for the case of a free electron (plane wave) with no 
scatterers present can readily be evaluated as 


~ l' tp 1 -JL£? _ 

s <Z\e 3 1 J e IA> 

- J _ 

= <4 1 C>j e 3 ju> = <4\e 1 u> 

r - e 3 \A > 


(2.14) 


Combining equations (2.1), (2.3) to (2.5) and 
(2.8) to (2.14) we have 


, -M t e 

a a p c 


p = e/sz^Ul^e * f/ttt A} 


where 

-/ _ -/t/y -pH ttH 

A = Pjj, i e [},e J e ia> 


( 2 . 15 . a ) 


(2. 15. b) 


The bar was introduced as a notation to emphasize 
the presence in our work of the configuration average. 
For the remainder of the work we shall omit this symbol 
with the understanding that configuration averages of A 
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and ft are implicit in all discussions and will be 
explicitly shown in all calculations. 
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III. EXPANSION PROCEDURE 
We cannot proceed to evaluate equation (2*15) 
directly. This would require the eigenfunction of H, 
i.e. the full solution to the N-body problem. Instead, 
A will be recast as an average of an exponential 
which, as Kubo has shown lends itself to a cluster- 
type cumulant expansion. 

We start with the matrix element, M, in (2.15) 


_ -it" -/3 H l*H y . 

Ms <Ale [>,e Je U> 


(3.1) 




= <4 le 


Define and by 

p id - pld 

G-Msl + fy (■*)= e*/*(-£*e 3 V f e 3 M 


2 ) 


(3.3) 


and the product G^M G^ (* ) by 

x ~ * 

GJ*)G A (y) s <2 £-jf e a C t ^ + 3 


(3.4) 


Note the identity 


aCb+c) = e ai e*r 


(3.5) 




^ 
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which is especially useful when Cb,cJ^ 0 • Upon 

substitution of equations (3# 3) to (3,5) in (3.2) 




o 4 rr 


M=<4le IP)] G, .(->'*) >2 > 


(3.6) 



Introduce the levelling operator L with the properties 


„ r ("(■*> , #=0,1 

L f t*> = i „ 

u l o , 


(3.8) 


so that 


r % 

6y^' ' 


Equation (3.7) can now be rewritten as 

-iu' £ o ( ,t/ , £ «,(/«» /«•<« - 

<?/" ' [ } t e ' f/" / Je <? ' W> 


(3.9) 


(3.10) 


Now define an average (signified by ( > ) for any 
product of functions of the three parameters ±it, 0 
by 

<f,«t)f,<PHi(-itf>S<jtie J f,t<V[),e f t (P)]e £(-‘VIJ> (3.11) 


so that M can be written as 

I Aw-*** 


n = <e t 


e; 


6 ' 


> 




N.V . > 


(3.12) 
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Finally introduce an ordering operator 0 such that 
even for non-commuting a and b, 
a t> 

e e * oe (3.13) 


and put M in the form 


> 


(3.14) 


In this notation Pj. , equation (2.16), becomes 

p JA = n> < 3 * 15 > 

where 1 is the identity operator. We have thus 
demonstrated that A is actually the average of an 
exponential. It will be useful for us to extend this 
average and define 



A clearly meets two of the requirements of being an 
average in the cumulant sense. Namely, Jl(l) = 1 
and the moment generating function is well defined. 
The remaining requirement, that the moments converge 
will be seen to be satisfied for all moments con- 
sidered herein. (A summary of the properties of cumu- 
lants used in this paper is contained in Appendix B. 
A more complete and rigorous treatment has been pre- 
sented by Kubo ). Therefore, we can write 
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= & %p 


{Jf(0e& CWnH/W+W*. i )} 


r e*/> {A‘(OLZ Tu (0, [?/‘V +%(/*) +%(■>*/])")} 


(3.17) 


where Ji means the cumulant corresponding to A • 
Further discussion regarding the cumulant A c in the 
specific context of this work can be found in Appendix 
C. The following examples demonstrate some properties 
of the simple (non-curoulant) average A* 

A (OL&ipSJM) = A 

_-/»<*“♦ V. i, ^ 

= Ciu < ': /, ' c L>i e je 10(3.1 a) 


A(PL2 3 (‘t>) = o 


(3.19) 


A (0L%(«>3,l-t )%,(-*)$ >p) % l‘t)3„ l-‘V ) 

-A(% i‘V3,M%(p)3, t- %, <• *>) 


-/jtf+W 

^ a/e R« Je M> (3.20) 


Equations (3.18) and (3.19) emphasize the fact that L 
operates on each class of g-funetions separately. This 



follows from equations (3.12) and (3.14). Equation 
(3.20) illustrates the properties of the g-functions 
as defined in equations (3.3) to (3.5). 

Substitution of equation (3.17) in equation (2.15) 

gives 


peZ-'U ( 3 . 21 ) 


This section started with an expression which 
contained the Hamiltonian for one electron interacting 
with N scattering centers. This has been recast in 
a form which involves instead a summation over n 
interactions, each involving one electron and n 
scattering centers. 
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XV. WEAK COUPLING 

In this section we briefly consider the first 
order of the cumulant expansion in the weak coupling 
limit. We will restrict our analysis to the time 
dependent terms and show that the Kubo formulation 
yields a relaxation time directly. 

Let us consider the time dependent portion l(tg) 
of equation (3.21) for n « 1 . 




(3 



g EM'dov) +A'(9/-‘V)] 


(4,1) 


Recalling the definition of the g- function in 
(3.3) and using the properties of cumulants we have 

5 ~ I * A(^tit)) - I (4.2) 


and 

= lA‘(3;(*)l 'f 


(4.3) 


where the asterisk denotes complex conjugate. 
We now evaluate the simple average 




xV 
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Jt(CjOt)) = %J<He e?p(-f.<t*e v,e^jfr,e * ]e IA> 


*<X\A>~‘ <4ie*/>(~S;<t*e * \f<? 1 )//> 


/Y 


-.£14 

*T“w „ a \ 


(4.4) 


Expanding the exponential in powers of V. which is 

J 

assumed small 


A (&<«>) 


, it fiU -£k it £% ~J 2 k _ 

z<2a>‘<XiJ-f.*e^V,e (t*e'v f e ‘)\....iX> 




+ (t <Xu »-’/<?/' la IV; I j'>l* - (4.5) 

Elimination of /l(Gj(it)) between equations 
(4.3) and (4.5) gives 


.^(V-,y) )=£ F Jd'iai 


(4.6) 


This integral is the same expression that enters the 
standard derivation of the transition probability in 
time dependent perturbation theory*" We can now intro- 
duce a* d&isity of states D(k). Since D(k) and 
K-AlVj'lTy can both be assumed to be slowly varying 
the integral simplifies considerably and we get 
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jC(fyty) t>(M \<X\VjiJ>l t (4.7) 

If all the interactions are equivalent, then 

(4.8) 

We substitute equations (4.2), and (4.8) and (2.l2.b) 
in (4.1) so that 

where 

r-'z^ypwl^ivi^l 1 (4.io) 

and 

f=A//I2 (4.11) 


Actually, the result just obtained is consider- 
ably more general than our earlier emphasis on mobil- 
ity might imply. The conclusions of equations (4.9) 
and (4.10) were arrived at independent of any consid- 
erations involving AH. It is, therefore, a quite 
general property of the linear response by a weakly 
coupled system that a steady state is approached via 
an exponential decay of the transient effects. This 
relaxation proceeds at a rate determined by T which 
thus constitutes the characteristic relaxation time 
for the system. 


t* r-» ■ i i*. , 




a- - •n 
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V. MULTIPLE HARD-CORE SCATTERING 
After a brief digression in Section IV where we 
considered the time dependence in the weak coupling 
limit we return our attention to the hard-core scatter- 
ing mobility. This section will consider some of the 
properties of the wave functions used to evaluate the 
matrix elements that appear in A* 

The most general matrix element to occur in the 
expansion of equation (3.21) is 
- itH(» .) ~{3H(n t ) ItH(ni) - 

<A\e [})£ le l-*> ) yi,+n x +yi t -)i (5.1) 


where 


(5.2) 


Vvxth the introduction of intermediate states and the 
judicious choice of formulae for the commutator we 
will see that much of the calculational work will 
require evaluation of matrices of the form 

-«H(n) _ 

<A\e \A> 


e <f t ie i 


(5.3) 


where 
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The hard-core scattering is represented by the 
potential 


f oo , 




L O , r>a (5.4) 

where a is the radius of the sphereical region from 

which the electron is restricted by the scatterer and 

for helium can be equated with the scattering length. 

-9 

For numerical considerations we will use a = 6.2x10 cm 
which is the value obtained by O'Malley in an extrapola- 
tion of effective-range scattering theory to zero ener- 


gy* 


1 


It is obvious from equation (5.3) that in order to 
proceed we must first be able to solve 


Hi ») if. ; r,,r ,, •■••?;)=£ (H ; fi •■■■ f») 


(5.5) 


An approximate solution of this equation, With H(n) 

defined by equations (5.2) and (5.4), was obtained by 

Coopersmith as part of the free energy study referred 

2 

to in Section I. The remainder of this section will 
present the Coopersmith method of solution as well as 
some of the simplifications which occur when it is used 
in a cumulant expansion. 

The approach underlying the solution is a replace- 
ment of equation (5.5) for hard-core interactions by an 
equivalent boundary value problem. We can alternately 
consider 
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(V x +-A') Tjf s ^ any f^><2 (5. 6. a) 

(S j f* "" ft) = 0 all f^<£a (5.6.b) 


This is a typical two region problem in wave mechanics# 
y satisfies the free particle equation everywhere in 
space except for n excluded spheres of radius a 
located at r, ♦ On and within the boundaries of these 
n spheres V must vanish# Because the ensemble aver- 
age makes the results independent of the particular 
choice of the n scatterers, we can for convenience, 
consider the set of scatterers as being labelled 1, 
2,#.. n. The Coopersmith solution then consists of a 
generalization of the simple scattering solution, i.e. 
a linear combination of a plane wave and n scattered 
waves 








h 

n 

t--! 


i-Att 


/>/«) 


jt. 


] 


(5.7) 


The A.(n) are determined from the n boundary con- 

J 

ditions that f be zero on the surface of each sphere. 

-4 ^ 

The choice of exp( IJt * ) for th© plane wave part is 

a matter of convenience and differs from exp(*«£*ft) 
only by a non-physical phase difference. Since every- 
thing is expressed solely in relative coordinates we 
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will uso the term coordinate to mean relative coord- 
inate unless specified otherwise. 

As indicated earlier, the significant features of 
low energy hard-core scattering are essentially deter- 
mined by the S-wave (lowest order partial wave) approx- 
imation. This suggests that we perform a multi-centered 
S-wave expression to determine the coefficients Aj(n). 
We first expand all terms in equation (5.7), except for 
(Afg,)' 1 expO A \r e , ), about r, , and equate the wave 
function to zero at . This process is repeated 

in turn about each of the remaining (n-1) scattering 

3 

sites and results in 

, A it* » / „ 

J. ' ***** +A;(») z + YL As») 1 lilh ( 5 . 8 ) 

A a. -/f cl J-i Jfc 


This is a set of n linear equations, one for each 
value of j. The prime on the summation indicates the 
exclusion of the term The solution of this 

system of equations in determinantal form is the ratio 
of two nxn determinants where the denominator has 
elements a gt given by 




Utr st 

6 sin fig 

IJia 

> 


S - t 


t 


(5. 9. a) 


st 


and the numerator has elements a 


given by 
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) t*/ 


u st - 


- e 


* A * ^r/ , ^ 




(5.9.b) 


In this form each Aj(n) is the sum of Ml 
terms. If we now consider 

<r c /e ir e ,> 

=/<** Yjliit.t. ) e'*%- («,; F, ... ?, ) ( 5 . 10 ) 


it is apparent that the combination Y'*Y' involves the 
sum of more than (nl) terms. We, therefore, will 
attempt to locate and eliminate as many of the non- 
contributing terms as possible before resorting to 
explicit evaluation of matrices. 

By inspecting the solution for Aj(n) it is seen 
that the determinants for A.(m), where m<n, are 
minors of those for A-(n). The minors are obtained 

J 

by covering up the rows and columns in the solution 
of A-(n) that contain the coordinates of the (n-m) 
scatterers that appear in H(n) but not H(m). Com- 
bining this with the ensemble average employed in A , 
it follows from the separability feature of ctroulants 
that all terms will vanish unless they contain all n 
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coordinates. This is equivalent to saying (in the 
standard language of diagrams) that unless these terms 
contain all n coordinates they will be unlinked. 

We also note that every element of the determi- 
nants which does contain scatterer coordinates has 
sin(ka) as a coefficient. This suggests a further 
expansion in powers of ka consistent with our mathe- 
matical formulation. The gaussian factor exp (-0k /2) 
in equation (5.1) limits the momentum states which con- 
tribute to k to the order of 0~^ or less. For 
ambient conditions this gives ka^ 10"** and for the 
conditions at which Levine and Sanders did their work 
(see Section I) ka^ 5x10 . With the expansion in 

ka equation (5.9) becomes, to lowest order 
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(Section VII ) that terms proportional to a 11 are 
also separable and the leading terms will be propor- 

Mil 

tional to a • 




■ 
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VI. SINGLE SCATTERING MOBILITY 
Before proceeding with the full solution for fA 
we restrict our attention in this section to the first 
order effect of hard-core scattering. 

Equation (3.21) for fi , with n » 1, is given by 



Our first task is to evaluate the various aver- 
ages. In first order cumulant averages are equiva- 
lent to simple averages, or 




( 6 . 2 ) 


To evaluate Ji we need the solution to the 
Schrddinger equation which for single scattering in 
the S-wave approximation is 


•,/v / j* . . r 

YX (Vty- (in) Ce -e 






nsr 


i 


(6.3) 


To second order in ka this becomes 


Vj (F. { Jj) = t> ri* - C - iM) 'l 


(6.4) 


*N 
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This leads to the diagonal matrix elements 

mM + Vi ) v 

si* A * a> 


. _ -i -«(^+VJ) 

= (wi/1) /-/V$ /r, A' <? < fi l ^ |ti-> 


z=(ev'n) '!■■■[ cfr e tfr c . e re > fjT %? r (r,;%)e'^ %,(£,;>)) 


-nra°<*2(mj‘‘d siJfzLj^jx e^e~^]e ~ 


(6.5) 


With the definition of Jl in equation (3.16) we 
obtain 

= '-^[(H't+DC (XT*) -t(K‘T- /jSOtfV ^(cc^-SiK ¥}] ( 6 . 6 ) 

£A‘(X<ri:-Rft+Mi<''*** L K)D(i</‘) -< WV-Oe^J} (6.7) 

I . i ' f 


where 


CM fly e.s ip- J | 


(Fresnel integrals) (6. 8* a) 


SW-ffy ^ J 

x 

(Dawson function) 

and the following dimensionless variables have 
introduced 


(6. 8*b) 

(6.9) 

been 


K‘=vf/J (6. 10. a) 

T= r//3 ,T=r //3 (6.10.b) 


•A ’ > 
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(6.10.c) 
(6. 10* d) 




A = a/ft'* 

R » a It pa ft 

The intervening steps for the above results are pre- 
sented in Appendix E. 

Expansion of the Fresnel integrals** shows that 
for K^T^> 1 the RHS of equation (6.6) approaches 
T/T where T = (2RAK)’ 1 * Thus we find that, just 
as in the case of weak coupling, a relaxation time has 
appeared as a direct consequence of the theory. 

If we limit ourselves to 2RA^ 1 (low density- 
high temperature; e.g. on$;afcfla.t and.'. 3O0nK) and assume 
that the gaussian factor restricts the contributing 
values of the momentum to K^l, then in these cir- 
cumstances Tg and the RHS of (6.7) is = 1. In 
these circumstances 







4 e 

3fV (zTTM 9) 1/4 


( 6 . 11 ) 


where m is the electron mass, 0 is the product of 
temperature and the Boltzmann constant and (T is 
the quantum-mechanical low energy hard sphere cross 
section^ •> . The normalization ’Z 

by Coopersmith to be given by Z=exp(-R). Equation 

(6.11) is exactly the same as the result obtained from 

A 

classical kinetic theory. 
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VXX. EXPANSION OF CUMULANTS 

fch 

Xn this section we will examine the n order 
cumulant expansion and consider the effects upon it 
arising from the levelling operation, cuiiuiant separ- 
ability and the thermodynamic limit. 

Each individual cumulant entering the expansion 
of equation (3.21) 

ol i | r ?/"< * vp > *$<-"> 1 H (7.D 

is necessarily of the form 


yit fi» f/i v 

A(0 LJX ft ypfi y-w) ) n, (7.2) 


>1. 




An essential feature of cumulants for the considera- 
tions in this section is that Ji c ($9 •) is expanded as a 
sum of products, each of the form • The 

prime on TT indicates that the product is subject to 
the restriction 2 = n • This is discussed in more 
detail in Appendix B. 

Let us first consider those terms where one of 
the tij =n and the other two are zero. The levelling 
operator will equate to zero any Jl which contains a 
repeated subscript. However, some terms will remain 


j 


< 




\ - 


'X 
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from a cumulant with repeated subscripts provided the 
repeated subscripts are distributed among the s 
such that there is no repetition of a subscript within 
a single J\ • (See Appendix C). As defined in equa- 
tion (3.16), sA has a normalization proportional to 

~ rf' • The ensemble average of introduces 

a factor proportional to S2“ >1/ + l • Thus each term of 

the n^“ order cumulant is proportional to SI ” n • As 
mentioned earlier, the dependence of Jl c on a partic- 
ular set of scatterers is eliminated when the config- 
uration average is taken. Thus, the n summations 
over the subscripts are just the number of ways of 
choosing the different subscripts out of N . This is 
N(N-l) ... (N-n+1) when all n subscripts differ, and 
N(N-l) • • • (N-m+l) when only m, (m<n), different subscripts 
appear. Now we combine the results for N and S2 above 
and go to the limit of large N and JT2 , N/O = p • 

For the case where all n subscripts differ we have 
(N/ft) n = f n • However, when there are only m differ- 
ent subscripts we have which in the limit of 

large N and _Q is vanishingly small. Such terms 
need not be considered further. 

It remains to consider the terms for which two 
(or all three) of the n j are non-zero. If one or 
more subscripts are repeated within a single set of 


V - - ■ - X 
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% subscripts there will be no contributions. This 
follows from the same reasoning as above and these 
terms will not be considered further. This leaves 
those terms in which there are only nU n different 
subscripts, but with no subscript repeated within any 
single nj • In this instance the same considerations 
used above show that there are terms proportional to 
/=> m . Since there are n scattering interactions 
involved, the non- separable (linked) term of lowest 
order in a must have at least n* factors of a • 

If n> m+1 these terms are of lower order in a than 
we consider here. If n=m+l, these terms are in gen- 
eral separable in the cumulant sense as is shown in 
Appendix H. The exceptions are 

and its complex conjugate with con- 
tributions^/ 3 ^ • is zero because 

it contains an odd function to be integrated over a 
sy^netric region. The former terms are considered 
further in Section VIII and Appendix I. The latter 
term is shown to vanish in Appendix J. As a conse- 
quence of the above we can rewrite equation (3.21) as 




~r 






r 1 A C 

* l Jt 
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+ it tyWftH + ty*/#-*/] )) } (?<3) 

Since most of the complexity of cumul^nt expansion 
is present for n=3 we have compiled all n=3 terms 
in figure 7.1. They are listed together with a summary 
of the highlights of the above discussion and refer- 
ences to the locations of additional relevant material. 

It has also been shown by Cooper smith that another 
simplification arises from the angular part of the in- 
tegration required to express the coordinate matrix 
elements in terms of the wave function, namely, from 

$ <Se sineYt ( 7 . 4 ) 

From equations (5.7) for Y and (5.11) for Aj we see 
that all coordinates enter multipli catively as 

either spherically symmetric scattered waves, 

(kr)"' 1 ' exp(±ikr), or plane waves, exp(±ik*r). Since 
the only angular dependence in Y f is in the plane 
waves, the integration of expression (7.4) transforms 
the plane wave part of Y* f into scattered waves. 
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J tM'V 

J o $ o St*# e. 


Ur -Ur 
e - e 

tJ?r 


Thus the exponentials in Y*V / can be written as a 
single exponential of the sums and differences of 
coordinates, exp(ik Z* P) ) . Considering the various 
ways of obtaining combination Z ± ^ leads to the 
conclusion that all such terms cancel unless they are 
±.ZT I £ , i.e. positive definite or negative definite. 
Since is real these will be expressable as sines 

and cosines of sums of coordinates. 

i 

There are two types of terms which occur in . 

They either arise from the product of a plane wave and 
a scattered wave or the product of two scattered waves. 
The latter terms must necessarily contain at least one 
coordinate in the exponential which has the opposite 
sign from the rest and will cancel those of the former 
terms in which the argument of the trigonometric func- 
tion is not positive definite. An explicit demonstra- 
tion is given in Appendix D for the case n=2 . 
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VIII. EVALUATION OF AVERAGES 
In this section we present the explicit results 
necessary for the evaluation of the exponential in 
equation (7.3). 




We will first consider the averages and then the summa- 
tions of the above expression. 

Let us introduce the notation Jl »(<*) to represent 
Jl (<*) ) and to represent that 

part of Jl(jj t J^) which is proportional to a m • In 
this notation 


A‘(L[t («> ( 8 . 2 ) 

From the arguments of Section V we know that 
m^n . Two values of m will actually be needed, 
namely m=n and m=n+l • The terms m=n alone are 
insufficient because (<* } a”) = O , as will be seen 
below. However, this does not make Jinfaja*) irrele- 
vant as it does contribute to the expansion of 






'N 


.A ' 
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Let us consider 


.x M 


u . . +E-V.) ~ni? x Hp x 

JU<V*(P„Sl, j S-fftmcZ^ * 'o.e *Je * M> (8.3) 




* 


= (n"<Xu>)' d\e ‘ t{ * + & V,) rt> (8i4) 


,» JL 


=si e f'-f%AAAe a>e ,M ‘ ifo ( 8 . 5 ) 


The contribution to from the coord- 

inate matrix element is 


tyir (-if a" rnfs A e *** 

0 •#! '!> *#•#,> i Mi/' 


1 


so that 

j cl *) s (" a rra ‘ ) 


h 


while the contribution to Ji H (if; a**') is 


and 


-i&f* 

^ * 


(8.6) 


(8.7) 


‘TirHfa”* , tt(n{MJl*e~ L *~ a*M (8-8>a) 

” 'at ht *he* 


47T6-/M ^ . f^7 (8.8.b) 
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so that 






Ji* 


All 

n*i - . a 


-/ </*[(«?''- on je 


(8.9) 


Equation (8.3) is as defined in equation 

(3.16). Equation (8.4) follows from equations (2.12) 
and (2.14). Equation (8.5) is the transformation from 
momentum to coordinate representation. Expressions 
(8.6) and (8.8) follow from the discussion of Section 
V. (Additional details may be found in Appendices D 
and F. (8.6) and (8.8) are extensions to general n 
of expressions (F.3) and (F.4) for n=3). The inte- 
grals in R.H.S. of equation (8.5), with integrands 
given by expression (8.6) and (8.8), are evaluated in 
Appendix G. 

In a similar manner we have 


H cXie^Ch + 




( 8 . 11 ) 
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^ t™" 1 * 4 '**) < s - 12 > 

The last average we have to evaluate is 

<*/**%, e'*** Vi (8.13) 

= <Aie itl >* V* **&» e l '*\ e* lM !t) ( s . u ) 


.* t* p^" rp'l. ^ 4 ^- . si a f lt i 


Alii 


= *^ 7 3 +fi'l(h(/»*')e 3 (8.15) 


Oj- >; 


The identity 1 for the commutator of an operator 
A with the density matrix 




(8.16) 


was used to go from equation (8.13) tp equation (8.14) 
where A is the time derivative of A. The transition 
from equation (8.14) to (8.15) is presented in detail 
in Appendix H. 

The remaining averages in expression (8.1) are 
obtained by complex conjugation, since 


\ - 


X 


(8.17) 
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A (JlWli (-■*>) *[JUW') ?/(/») ]* (8.18) 

We note that the cumulant separability of 

for n>l follows from equations (8.8) 
and (8.11) which show that 

=[A(«i *)]" (8. 19. a) 

We also note that 


sAfltti*) +sfl c t - o 


(8. 19.b) 




(8. 19. c) 


It now remains to perform the various summations 
indicated in (8.1). As discussed in Section VII the 
n-fold summation over j from 1 to N, gives 
N(N-l). . . (N-n+1) which goes to N n in the limit of 
large N . Thus, performing these summations has the 
effect of replacing Jl“ n by f>* in (8.8), (8.11), 
(8.14), and (8.19.c). 

The final procedure remaining is the summation 
over n of the cumulant expansion. 


= it WJ 2T (-// (f M! 

**! 



f, 

?n * 


( 8 . 20 ) 
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( 8 . 21 ) 

( 8 . 22 ) 


££ 

pltf/l) 


similarly 



In writing equation (8.20), the standard expression 
for the cumulant expansion (see App. B) has been modi- 
fied by separating out one A to demonstrate explicitly 
the proportionality to a*"' 1 ' and inserting P to 
indicate that we must include all the permutations of 
the remaining A ' s that satisfy this proportionality. 
Equation (8.21) is a rearrangement of the summation 
which contains explicitly the combinatorial factor 
arising from P and also makes use of equation (8.19). 
Equation (8.22) follows from equation (8.21) as a con- 
sequence of the properties of the Stirling number of 
the second kind (see Appendix H). Equation (8.23) is 
a straighforward summation of equation (8.22) with the 
Jl's eliminated by using equations (8.7) and (8.9). 



r . <r* x m-t) 

[+ 'Xit-Xfc -2 -ae ] (8. 23. a) 


< 



The remaining cumulant is 


f Jtty k>>) - tl$ (nr>) -M/I'>)A(%II»)] ( 8 . 24 ) 

The second term of R.H.S. is smaller than the first by 
a factor of since each vfl is^A" 3 - and there is 

only a single summation ^ N . Therefore t the second 
term is negligible in the limit of large N and A . 

In the preceding sections we saw that the cumulant 
expansion is equivalent to an expansion in y° and the 
S-wave approximation with its related expansion in ka 
introduced a further expansion into each cumulant. 

Our approximation scheme now consists of locating the 
leading term in a for each power of f 3 and summing 
these over all powers of p . In this sense we are 
solving the problem tu : ’all orders. 1 ' 
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IX. RESULTS AND DISCUSSION 
As a consequence of the work in Sections V,' VII 
and VIII, equation (3.21) for the mobility can be ex- 
pressed in terms of cumulant averages by 

ft = 4 ^ 2 '' U [J,'e 


x e 




(9.1) 


Upon replacement of the cumulant averages by the explicit 
expressions obtained in Section VIII, we have 

^Swn k e e l<* Te ■ (9 ‘ 2) 


where 


(9.3) 




? -WJltlV*'"') 

= A R { ( r'+ r‘) DM -(a'“+cc‘) p(a.) +(Q t -2<z‘-<i'/ , )i>((Z) 
„ +J% : [aRK '‘G:‘e^ - 0 - o/ ')e^ ]) 


(9.4) 


and f p is the time independent part of the second 
term on the RHS of equation (8.1) 
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.jv. ■£ r' a u* 

'* Rf\e JTW' 

Similarly 


e a 


f~f^f x 




(9.5) 


(9.6) 


= 2 /? /? { r<? t [ S ( <?. )* C (Q.fiZ) ] 

**?[$(* m) -C(K^W)]- ti / '[S(a</Vr')-C(Q t JW)] 

- Q. '‘[(sihRT - CCS RT)C(Q.Jfir) + (St'ngT +osgT) S(0-f%r)] 

-JVn(s‘»ToS- cctTO;)] (9>7) 


$ '.a _j 

F x S A (a. 71)'* ft /) e J c </* (fi+A) * 6 “* ( 9 . 8 ) 

and 

Qi s K' + sR (9.9) 

Equation (9.4) follows from (8.22) and (8.9) while (9.7) 
follows from (8.22) and (8.7). 

It is interesting to note that for all physically 
reasonable values of the parameters f^ and E^^F^. 

Thus no eumulant with repeated subscripts contributes to 
the results. 


43 


When we expand C(x) and S(x) asymptotically we find 
that for Q 2 -€ T 

FSr U. R A ( R + K x ) T (9.10) 

In general there is a large region in which this is sat- 
isfied. Typical values for a mobility experiment are 
lO" 1 ^ 2RA(R + K 2 )^< IQ 2 and 10 7 < T £ < 10 10 . Thus we 
again (as in Sections IV and VI) find that after an in- 
terval small compared to the duration of the experiment, 
the transient effects die out and the system approaches 
a steady state exponentially. We therefore make the 
identification off the relaxation time as 




(9.11) 


We note that where K»R ( e .g. one atm. and 300° K) equa- 
tion (9.11) reduces to the first order result 
T = (2RAK)” 1 . As before, this relaxation time result 
was obtained independent of the specific form of £H. 

Figure 9.1 shows exp(f) as a function of K for 
selected values of R. In each instance the function 
looks like a step function with the step located vary 
near K=0 and the step height dependent on R. For R^l 
(e.g. one atm. and 300° K), exp(f)=l, as was the case in 
the first order calculations (see Section IV). The inte- 
grations indicated in equation (9.2) were performed nu- 
merically using a Simpson 9 s rut© for the K integration 


i 
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and a Gauss-Laguerre procedure for the T integration. 

In the former it was necessary to assume that the gauss- 
ian factor exp(-K /2) provided an adequate cut-off, as 
exp(f) becomes strongly divergent for large K. This is 
a reasonable restriction as the inclusion of large K 
is equivalent to considering contributions from high 
energy electrons which is physically unrealistic and con- 
trary to our mathematical restriction KAtf 1 . 

Figure 9.2 shows I* vs. p for a temperature of 
3.96° K. The circles are , the experimental results ob- 
tained by Levine and Sanders. The dashed straight line 
is the semi-classical theory given by equation (6.11). 

The broken line is the first order result obtained by 
numerical integration of equation (6.1) with the Jl 1 s 
expressed by (6.6) and (6.8). The solid line is the 
full result obtained by numerical integration of (9.2). 
Since the integral of (9.2) is a function of the dimen- 
sionless parameter R = 2Ff>ap, it follows that the solu- 
tion at any other temperature (i.e. p) is a simple 
translation (in/ 3 ) of the above result. 

At lower densities the experimental mobility 
approaches that predicted by the theoretical calculations 
for free energy. The high density region is understood 
in terms of the formation of a correlated, or "bubble", 
state. However, the transition region in which there is 
a large drop in f4 fo?: a fairly small increase in f has 




not been explained* Calculations based on a transition 
from free electron motion to correlated motion predict 
a transition that is much too sharp* Our calculation 
indicates a decrease in /U. that is qualitatively cor- 
rect (i.e. the theoretical and experimental curves 
have similar slopes) but the result is quantitatively 
wrong (i.e. it occurs at too low a density). We do not 
see any modification of the theory or improvement on the 
hard- sphere model that could cause the theory to coincide 
quantitatively with the experiments. 

It is rather interesting to find that our theory, 

which considers only free states, yields a mobility which 

becomes! vanishingly small. The severe drop in only 

appears when all higher order (multiple scattering) terms 

of the cumulant expansion are included and comes from 

exp(f) (see figure 9.1) which in turn comes from 

(AH,e"^J • In the formulation of the theory this is 

the interaction of the external perturbing force with 

the equilibrium (dH=0) density matrix. However, the 

physical process causing this rapid change in f* is not 

clear. One possibility is that our result is somehow re- 

2 

lated to the Anderson transition whereby conduction goes 
to zero in certain random systems. Another possibility 
might be the appearance of some ordering in the system. 

t 

Both of these suggestions are highly speculative, but do 
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indicate directions in which, this study might be pursued 
further. 
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X. CONCLUSION 

We have applied cumulant theory to Kubo’s formula- 
tion of irreversible quantum statistical mechanics. In 
the weak coupling limit we saw (Section IV) that in low- 
est order the approach to steady state is always charac- 
terized by a relaxation time. Similarly, for hard core 

!> 

interactions we found that with the use of Coopers£<ith 
wave functions we derived a characteristic relaxation 
time independent of the specific nature of the inter- 
action, /6H. For the specific application to mobility, 

</\H = -zeE, we found that the general theory yielded 
the well known results of semi-classical theory at high 
temperature- low density. For the region of lower temper- 
ature and higher density the theory was shown to give 
qualitative agreement with the anomalous results of ex- 
periment. However, the quantitative disagreement indi- 
cated a need for further work on this problem. 
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APPENDIX A 

KUBO’S LINEARIZED ADIABATIC SOLUTION 
OF THE LIOUVILLE EQUATION 1 ’ 

The density operator y® T (in the notation of 
Section II) obeys the Liouville Equation 

= ,n] (A.D 

This is expanded as 


= V] * [»> + O H,c] * [ (A . 2 ; 

The equilibrium density matrix p satisfies 

4 ir - C H > f] 

and in a linear approximation we can neglect 
Thus 


+[H,Ur ] 


(A.3) 


We now solve for Af 
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-WAf * * jfcAf + Af>H -l AH if 3 ] 


- i'fW ,'ftf -lYM ,'YW 

e if (e Af>e )e 


<' [AHiP] 


itH -itH 

e Afe 


e [ AH>;]e 


where t is that time at which = f ^ Afft.) 
This gives 



e [AH,r1 


e 


which with the variable changes (t-t f )-*t and 
(t-t )-» t„ , becomes 

O E 7 



-itH HH 

e e 




\ ~ 


(A.4) 

o 


(A.5) 


(A. 6) 
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APPENDIX B 

CUMULANTS - GENERAL THEORY 1 
This appendix summarizes the properties of cumu- 
lant theory that are used in this work. Cumulants are 
also sometimes referred to as semi- invariants or (in 
complex form) characteristics. 

The defining equation for the cumulants of x is 

JL‘(e**-l)- (b.i) 

where J\ is a normalized average and Jf is the cumu- 
lant average corresponding to the simple (non-cumulant) 
average Jl . Both sides of the equation (B.I) can be 
expanded and coefficients of like powers of f equated. 
This leads directly to 


A‘0)=A(D--l 

(B.2) 

A ‘M -- A M 

(B.3) 

A' (■*')* A(x')-A‘M 

(B.4) 


(B.5) 
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S!‘(zV =S!tx'')-[vJt(x , )JII*) + 3jl'h'l]*iiJI(x‘)Jl‘(x)-(J) v <*) (B.6) 

and so on. 

Straightforward extension to multiple variable gives 

Jl c (&**-!) xj^JUe Tif ^) (b.7) 

which yields the following relationships: 




(B.8) 


JV (v,7iXi) = JK-x.tiXi ) - +J1 (y,),/Ux i x>) +Jl(Xi)Jl()(,x t )) 


+ aJUxOjWy^J? c*i) 


(B.9) 

— * 

and so on. 

A general expression for ** "’*>i ) where no 

variable is present to a power greater than one is 


M*" IfOlY. fr J\ (jf *%) 

f*?) Jr. r** 


(B. 10) 


2 

ft! * 


The summation in (B.10) means the sum over all possible 
arrangements of the set of n variables in q subsets. 

Much of the utility of cumulant theory in many-body 
studies comes from its ” separability" property which can 
be stated as follows: A cumulant will be zero if one 

(or more) of its variables is statistically independent 
of the other variables. 
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As an example of separability we consider equation 
(B*9) with Xp separable* Separability of x^, implies 

\/l y % ) * >/} (V|) v/?(V») (Ball* &) 

JU*') (B.U.b) 

iA (v.VtXi) * (B.lJ,.c) 

Using equations (Ball) in equation (B*9) gives 

i/2 x» x») z ,/i(Yi)>/l(y* Vi) * [\A (y>)JUy\ v>j # ( , X\)\/l(yi)\/H)t)) +\fl (ty)\A(y>)*A(fa) ] 

+* a (y>)\A fa*) (xi) 
s 0 


(B* 12) 
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APPENDIX C 

CUMULATES' - IN THE PRESENT CONTEXT 
The four step prescription for disentangling equa- 
tion (3.17) is; expand the exponential, order the re- 
sulting products, expand the cumulants in terms of the 
corresponding simple averages, and lastly, apply the 
levelling operator. 

A specific example of this procedure is 


X 9 






(c.l) 


where 

$s 9 Hi) +9(p) + 


Dfc&ajase; 
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A typical expression that would require both order- 
ing and levelling might be found in the third term of RHS 
of equation (C.l) where the ordering would consist of 

A‘(0l%U»1,W3M) ' A‘(L JJftl ?; M) (C.2) 

The RHS of equation (C.2) can be expanded in terms of 
simple averages by 

= A (L -[A(L tw)JHL 9>>) * VM(L *,<*/?,</»)] 

+ (0.3) 

Finally, application of L to RHS of (C.3) gives 

A‘(0L r -s[A(%U»)A(t m) (c.4) 

In general, the operator L must be applied after 
the cumulant has been expanded. In the above example 
this means that L can be applied to RHS of equation 
(C.3) but not to RHS of equation (C.2). Otherwise the 
surviving terms in RHS of equation (C.4) would be lost. 

In the present work, however, one could safely commute 
the cumulant expansion and levelling procedures because 
the surviving terms, as in RHS of equation (C.4), turn 
out to be negligibly small. This last point is discussed 
in Section VI. 




xV 
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APFENDIX D 

TRIGONOMETRIC FUNCTIONS OF COORDINATES 
We noted in Section VII that the arguments of the 
trigonometric functions entering the coordinate matrix 
elements are all positive definite# We consider this 
matter further in this appendix which contains an expli- 
cit demonstration of the feature for n=s2» 

Before proceeding let us introduce some notation and 
collect some relations found elsewhere but needed for this 
presentation 
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e sen 9 e 


<‘Ar -<A» 

q - e 

i 


Q.< ~ ^ it 


(0.5) 


We will also use the following schematic notation for 
locating terms in the determinants. A circle at the 
location of an element indicates that the element contri- 
butes to the result when the determinant is expanded. 

The angular integration of 


%*(£ ; r,,P,) e ^ ) 


(D.6) 


is 


f do sen 0 [e *' +/), (*) *A X ** +A\( % )dgu J (D. 7 ) 


Let us consider the particular cross product 

J <to ti*e e f 


( 0 . 8 ) 


and look for those terms containing indices 1 and 2 and 
proportional to a' and a (see Section VII for fur- 
ther discussion of this point). The terms we want are 


lame* 


x\ 
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"X* S"* e I 771* 77f i e '' 


r K , -tJ'Y,* 

~-\JesUee /(-*„<? * (/+*, JM, a *) ja ( 


»/« 


*v 


-l(4&) + (&*“&€! ] (D*9) 


let us look at the exponentials contained in R.H.S. 






(D. 3,0) 


CLt\a n a e ,r=* € 


t A (~ ha ^ ha * h"/ ) 


(d. 11) 


a iA(r t , + 2 Kr + Yeu) 

a ei a u &*•! £ 


(D. 12) 




lAl-K' + alrii + fc.,) 


(D. 13) 


The coordinates in (D.10) and (D.12) are positive 
definite while (D.ll) and (D.13) are not. We will locate 
the complex conjugate to (D.10) and show how it combines 
and also the counterpart to (D.ll) to show that it cancels 
Let us look for a* 2 Ctn • The only other 

source for a €t is from the scattered wave A* l*) <2 #a 
or 

a it (A' y e ', # 

J e si »0 e A 3 


(D. 14) 
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which gives 


, * o I 

iSlLL Z 9LslL l I ¥■ 

aa r oi " 

from which we can pick out 

( t d 6.) & 4 i & 1 (l + + ) 

The other cross-product of 

\jt> /)*M e‘ x /I, (*)#,., 


which gives 



from which we can pick out 


(.•>«)' V, a„a t ,, - x a,’ <?„ a r , 


(D. 15) 
(D. 16) 

this scattered wave is 

(D. 17 ) 

\.'\ 

TJ ) *<•! (D. 18) 

(D. 19) 
of R.H.S of (D.9) with 


Combining the first term 
(D. 16) and (D.19) we find that the exponentials with 
ordinates of mixed signs cancel, and we are left with 


2 Ci U ( hi. * Ki'* ^f 1 /) ” & * K \x t ?*•/)] 


(D. 20) 


59 


APPENDIX E 

INTEGRALS, FIRST ORDER 

Equations (7.6) and (7.7) are obtained as follows: 


Jt t s £ [A‘ ( S'"' 1 ) + A Y f { !-<*>)] 

*£[Jl<9"'0 *<a (%<■“>) ] 


r< 


= £[ I >A(6j(-‘0) - I ] 


_ .•‘*(r*U -ae iff 
" ,3 D,e 3 ]e 3 U> 


„ *£?' P i 

+ <Me*[>,e Je i1>-zZUJ 

»> * e%^U a / 


Similarly 


/ ~h 

-Z r fa £%!<*(> (tie fl( * ' jj> - /] 






»iai*sa 


vV 
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Thus we need to evaluate 


°Ua { 


m <ur K ' ,v ‘\x> 

<<? iK> Oiie ^ \ i>/> <?,'</> 


> 


% A t i+t At* m £$*(%* ^ 9 *) 

=faff)M/A/^'e <V5i42 / f>/ 


where 


<>*/ « i»j<> 








.. p"'^V, -it-< -•/•fv . _ . . 

5 (iffj i e - 1« ^prj- j 5 - ie ' VM<6%7 J 


(inf'f/A e * e 


-«^ a - t ' 4 '(K-fh') 


**^(%/* % j) 4"'S Cc* tort 


-iU' 


-J&±£ a r 

M' A e a " 


-<2 


a«J 


I 1 /*/} 

+W'l '* 1 


4 
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Eliminating tha matrix element in ^ we have 


~nau>e 





or 


<* M =0.e 3 {<-fiZ> 


uo„ 

(**?”* 



With the substitution of °<^ for the momentum diagonal 
matrix elements of Jh and Jl/% above, equations (7.6) 
and (7.7) are obtained directly. 
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APPENDIX F 

COMBINATORIAL FACTORS 

In this appendix we show how to locate the parts of 
the expanded determinants which enter into equations 
(8.2). This involves the explicit functional form of the 
integrand as well as the combinatorial factor associated 
with the number of ways of obtaining equivalent terns. 

It will be convenient to introduce the following 
simple diagrammatic notation for the terms to be consid- 
ered. <X st > stt is represented by a straight line 
labelled s and t at the two ends. Repeated subscripts 
are attached their common point. a gg is represented 

by a dashed loop attached to the point labelled s • In 

2 

this notation the term proportional to a in equation 
(D.10) is shown in Fig. F.l.a. Those proportional to 
a in equations (D.10) and (D.12) are shown in Fig. 

F.2.a. Complete consideration of n=2 would have turned 
up the additional terms indicated by the diagrams in Fig- 
ures F.l.b and F.2.b. 

The equivalence between the determinant notation of 
Appendix D and the diagram notation above is shown expli- 
citly for n=3 in figures F.3 to F.10. The determinant 
shown is the one resulting from the angular integration 


sssrace 




<!.■*?** ' 
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J [^<S> sn\$ e iX **'* A 3 (3) a tt (F.l) 

The general structure of as the ratio of 

determinants leads to the following conclusions. There 
are nl equivalent terms and having all n scat- 

tering indices. This corresponds to the nl ways of 
permuting the n indices among the n vertices exclud- 
ing the fixed end points, e and e» (see figures 
F.3.b and F.4.b). There are two classes of terms rs cc 
having all n scattering indices. There are the terms 
identical to those above except for the addition of a 
dashed loop. As there are n vertices available for 
the loop, there are n(n I) such terms (see figures F.3.a 
and F.4.a). Finally there are the terms which consist of 
a line of n* vertices between e and e» with an 
attached loop containing (n-n* ) vertices aside from the 
one shared with the line. For each n' such that 
1^ n'4 n-1, there are n'(nl) equivalent diagrams. This 
can be seen as nl was of arranging the vertices, given 
a specific shared vertex for the loop and line, and n' 
choices for the shared vertex (see figures F.5.b to 
F.lO.b). 

The diagrams of figures (F.3.b) and (F.4.b) repre- 
sent equivalent integrals, as do those of (F.3.c) and 
(F.4. c) , (F. 5) to (F.8) and (F.9) and (F.10). With the 
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aid of these figures, and recalling the definitions of 
equation (5. 11. a) we see that 




Ui * TSiktfXf.) 


contributes 


~«A 


-«na'3‘j M4e 2 

0 N? I 'll rn If}#* 


to A and 




-una* 3 l[ e 3 4 + 

J * L to,tun»n** 


(P.2) 


(F.3) 


•f 2 


linjji ft,-* IhJJkd 

ht Ijt Kc' 


+ 


s<‘v> A ( 1*1 + Oil !j± t \ljJLhill 
ir« h* h» r„ 


] 


(F.4) 


to Jl, (* ; 
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APPENDIX G 

INTEGRALS - ALL INDICES DISTINCT; n th ORDER 
This appendix contains the integrations necessary to 
go from equation (8.5) to (8,7) and (8.9). The presenta- 
tion is based on the work done by Coopersmith in his study 
of the free energy.' 1 ' There the essential calculation was 




e 


rw 




-fiHl*) 

e l&> 


where the trace was evaluated in coordinate representa- 
tion. In the course of that analysis it was necessary to 
utilize the Fourier Integral Theorem. In our notation 
this is equivalent to 



The bracketed expression appears as the basic building 
block for our present work. Therefore, we are able to 
adapt his solutions provided we are careful to select it 
in a form in which the final step (i.e./V/ ) has not been 
taken. 
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Let us proceed to evaluate the integral obtained 
from equations (8.5) and (8.6). 

I , ,< c 

( " 4 ‘ ’ 

where S = ( *«» + **»*• * ' " + ' 


First we consider 


-W 


*~ s ' s 


V x 


(JT)' 1 1 

a 


- s 
T7t 


- iA’l^e “ ) 

».•■* e . $ e 


This gives 

r, = ( ^nr ^ ^ T* ^ ^ ^ - w 

B=2 e ~ r (-11 f*)" « * 


X 

T7T 


As there are 
variables of integration we 
and redefine s 


n+1 relative coordinates and n+2 

can transform to new variables 


*-£»■/ 


-_C 

-/ *!+«)/» /*! r- J)S r fi a< s iVi-^K 

(?, - a. v? J;J # ^> se 


The method of solution is as follows: It is clear 
that the effect of the n + l integration could be ex- 
pressed equally well by using a linear differential oper- 

at or i such that 






A 


"A 
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= 0H Wi'iiA */a) 

It is therefore sufficient to find the effect of perform- 
ing this operation. An equivalent expression for 




k exp(-ttk /2) was seen above to be 


J £ 


-i*d l 


- s 


a ~ 


[irr) 


"l 


fitpdS s s e 


so' that we can let Ms s-r, and write 


Rh* B M 4 ml f 0 dr, [y* s e scn^r, if sc*4m. 


- s 


= i ^ K' 1 fj 1 [Stif'S* s e * (t s < « * (vf'AJt 




= & h fi” (**>"' {* f 0 s 6 *'* {“>s[* s -< A +v k ]- £ * i [4s + (l''J)lr]'l 


- s 


- s 


s B* {i& K’ 1 i? I s e * ‘ * (4 V siH * s ~ k^Js) 


zB 




n 2 




„ ,/TT „ -xj\, -±<* X S) v 

^ to I. it l -* e ( fer— ) 


= ft 


(21 T) 




.kit 


... . -A'ui-y) -*1? 

B ™ -<# ■ M~ f yx { 4 e e 

2 k-*a L 4 


rS^:~rto. .tosto= ^r-SMMese^^ 


n\ 
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r» 


Contemplating the initial expression defining •£ 
we note that 


./•/ 

*- itt) 


i'Je * =P H U] 


Thus, 


ft = bJ^l' p P, U\ it-*') 


However, we also note that 
/TP 

jf dX X, Pjd) it-x'l =r <?„U , it) 


where Q n is also a polynomial with n terms. But 

Q , must equal P : therefore P., can contain only one 
^n-l n n 

non-zero term and from explicit consideration of P 2 (not 
shown here; see reference (1)) we see that this is the 
highest power of (it), namely a n (it)^"^* • This gives 
the recursion relation 


a 


h ~ 


2 


so th£t we have 

-itA 1 

f% = a j •» yu 


( - a n f a c t) 


M 


In a similar manner we now evaluate the integral 
arising from equations (8.5) and (8. 8. a) 
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„ , v/ . .iftH’K') „ ,iU“ 


.(U‘ 

j/Jj'e = 


m y 

<^ % (* !Qj a ztt . (zv) . f r ,'*-jL. -f’t/’ 

Ts ~Z/i)* e T< 7*1*1* JI* 1 


- S 

ITT 


-£A>(\r t -*i>) jpT* 

S » ' (.wfrsT I"'!/, ^ ' «>„••• t*.fc7 e 


l'<v*‘ 


l II iu <) i\ 

:'f)-V)' K « <s a (- ff W MHH***]#- 1 ) 


C, S = <‘' f > 


S,. c„ -T *”/./ SV? « '" ...A 


-> s 
37 ? 


C„ C l ~Ss e ** s^W* 


<, <*j* //" e “*■ [ d, * a 


.«£? , w^/ 

c»u'*)**l e L ^ e 


X 

a 


fl & . ,*+i , 

r £„ (Li)''l ch f A J e 


-X at -y) 


ft? 


- c* 


M "4 (it ) 

(Yf _ XV £ 

" »! 


zssr. 







where now s is defined by 
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- D* - i 1 ( A f £ a** « € 
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APPENDIX H 

SINGLE REPEATED INDEX; n>2 
We show in this appendix that where there are n 
scattering interactions and only m«n~l different scat- 
tering centers the corresponding <A is separable to 
order a n . 

Let us consider a set {ml of m different indices 
divided into three subsets f m^j , I m 2 ? and fm^J such 
that { ntj) and l have one index in common. Thus f m^i 
contains the indices 1, 2. • .V . , .m^ ; m 2 contains m^+1, 
mj+2, •••• m^+rr^; m^ contains m^+n^+I, m^*fm 2 + 2 . . . .m, 0 . 
This particular choice is convenient but not necessary 
for the results to be obtained. Specifically, we have 
to consider that part of I which has all n inter- 
actions represented and is proportional to a n . 1 is 

defined as 



s^T 0 


lo is represented diagramatieally in figure H.l.a. The 
topographical factor associated with how many ways we can 


obtain this integral for a given m is 


m nun* 

2 (rtf, -Of 2 1 yyi t !(w>-i)!i: 

‘ c ~~ V OTi ,)->«„ 






can be rewritten as 


lip Rtf'S e 


where 


V s <Kie 




l ?„•> 


a 1 t> -tA'(K**t K*t*) , i 

h’skjf&'drNfl'f e VI*XWU 


W*i 


x e 

- /£• (£*+»“ ^»*v) . 

x fi VKWU^ 


lr x ‘dr HP 1 ) e 


a-t 




m (^ '(fm+r 4x > . x 

«<? VO'tf^V, 




a/*** ^ . 

- i A h »+* ) 


w/»> /*•*). 


X<? 


X £ 


*»fi, XKV 
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whore l m^{ contains the indices mj+mn +2,....w,0» 

In the last step we introduced a Fourier transform 
(•'opened the ring") so that is now indicated diagram- 

matieally in figure H*l.b. This figure also shows the 
transformation to relative coordinates which gives 


i . tfc i# i<V‘> i,<4T) h, a r> 

‘'hM WtjH ., lXs 


ZJW)Z IJ 0, Hjl, A e‘ H '^ *V(fis 

T U 1" I" T" f ffr Jt ^ ^ ^ 

J S J4J)= jp* £ 




IS, 


The coordinate matrix elements are given by equation 
(7.6) which is this notation is V(it,fs})_ . , so that 

6} B 


r H/*q S> W/Vs»-0/ f" C& n. ct 

k ~ iTi s * m '(* "*U<') Sl J o * t* > ^ > w ' v« # 


3/i 


where K'= A'-Jt" 


This multiple integral is evaluated in Appendix G 

(as part of the evaluation of R ) with the result 

n 

SS -y/7 * £ 
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Similarly 


vXWHftK* fit i 0 ^ 


-JS &41 


=-Z 




17 v ' l " A U-i)! ($,-■>/ '0! 


e 


_ ~V-rr(srx)l('-xnci) 
jpp— — 


1 X Zi£L!£L& 

I P \6 P "T* _ . jTTM’ 4- /f. . .i it j / 


5,-V 


-//>*’ 


where k « (k-k") 


The evaluation of 1 is completely analogous to 

s 3 

l g with s^ • replaced by Sg and (it) replaced by 
(-it). When the above integrals are collected together 
we have 


J zi-zrrait/ji)* 1 '’ 1 (-ina fl/fL)*' (-2TTa it ) W,W 


_££ 


E -4V-^V. 

— Jt** ,'jL> fjr'ILzJL )g^=* (l~ e - 

n)* A'-A * }** I / e l Jr'--* 1 / 


Comparison of the coefficients on the upper line 
with equations (8,7) and (8.11) shows that the. simply 
linked indices factorize and are therefore cumulantly 
separable. The lower line is the contribution from the 
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repeated index and the integral is the reclosing of the 
link at V . 


77 


APPENDIX I 

SINGLE REPEATED SCATTERING INTERACTION ** P*? 

In this appendix we look at that part of 

2 

which is proportional to a and contains the dual 
effect of Vj . This immediately restricts us to the 
product of functions of (it) and 0, each proportional 
to one power of a . 

For this appendix we define A = A 


(tie'*** * V/ \he fl 1 1 * Vi) )e^\X> 


(i.D 
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-*''*'*. r‘i, r fi j . J'<*- 


0,f%e r( *' f> p4e 


(Jrr/% I- cr * 1 - "T 6 




(x.7) 


*f t/h-i+x ) "* 


(1.8) 


L f 


g a 

(/*+/*>** 
-Alt 


(1.9) 


4l£ 


fffl, (tj f e *_ ^ f l j c g 3 1 

*<**)* lit uS'r* n /f (^/>) j/x J 


(I. 10) 


* o» -/»>’* ^ (7^7P‘ J 


(X. 11) 


In equation (1.2) we have introduced an identity 
for the commutator.*' Equation (1.4) is obtained by 


taking one power of a from V(it) 


kk. 


and one power 


of a from the product V(p). - V(-*). w V(a). . . 

The last term in (1.5) vanishes because of k^ z which 
makes the integrand an odd function over the symmetric 
region of integration. The order of integration is 
changed in going from (1.6) to (1.7). (1.9) follows 

from (1.8) with the change of variables x-y=p, y=q . 

The complex integration in (1.9) is taken along the real 
axis from (-q) to zero and along the imaginary axis from 
zero to (it). A further interchange of the order of in- 
tegrations gives (I. 11). 
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APPENDIX J 

SINGLE REPEATED SCATTERING INTERACTION; 
NON-CONTRIBUTING 

We use the same technique as in Appendix I and 

2 

select the term proportional to a which contains both 
the (it) and (-it) aspects of the interaction. 




*[>,e ]e u> 

i 

-- < V(‘%. (A; ’ ) V(-«U‘A 


The above is clearly seen to vanish as soon as we 


recall (see Appendix I) that, in order a, v'(a) 


kk' 


is 


an even function of k and k' . This makes the integral 
over d£' vanish as a result of the oddness of its inte- 
grand. 
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APPENDIX K 

STIRLING NUMBER 1 OF THE SECOND KIND 
The multinomial coefficient M is the number of 

q 

ways of partitioning different objects into m^ subsets 
e'ach ^containing k objects 

mAifrw 9 Ifo 

The Stirling number of the Second Kind S^^ is 

q 

defined as 


— 


M 


zr^f 


«- 


It satisfies the summation relation 


£ s 4J =/ 

/=* ip 


Thus from equation (8.21) we have 
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APPENDIX L 
CUMULANT SUMMATION 

Equation (8.23) is obtained from (8.22) as follows: 
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None (~p8 3 ) 


Location of discussion 
and calculation 

Section VI, VII, VIII 
Appendices E and F 

Section VI 


Section VI 


Section VI 
Appendix I 

Section VI 
Appendix I 

Section VI 


Figure 7. 1 - A listing of all n = 3 terms arising from the expansion of p by cumulants. 





Figure F. 1 (a) Diagram corresponding 
to equation (D 12). (b) Equivalent 
diagram for n = 2, noi considered in 
Appendix D. 
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Figure F.2 - (a) Diagrams corresponding to 
equations (D. 11) and (D. 13). (b) Equiva- 
lent diagrams for n = 2, not considered 
in Appendix D. 
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Figure F.3 - (a) Schematic notation whereby the elements of the determinant 
which contibute to equation (E.4) are indicated by encirclement, (b) Equi- 
valent diagram notation for term ~a . (c) Equivalent diagram notation for 
term ~a . 
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Figure F.4 - (a) Schematic notation whereby the elements of the determinant 
which contibute to equation (E.4) are indicated by encirclement, (b) Equi- 
valent diagram notation for term ~a . (c) Equivalent diagram notation for 
term ~a^ 





Figure F.5 - (a) Schematic notation whereby the elements of the determinant 
which contibute to equation (E.4) are indicated by encirclement, (b) Equi- 
valent diagram notation for term ~a . 
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Figure F.7 - (a) Schematic notation whereby the elements of the determinant 
which contibute to equation (E,4) are indicated by encirclement, (b) Equi- 
valent diagram notation for term -dr. 
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Figure F. 8 - (a) Schematic notation whereby the elements of the determinant 
which contibute to equation (E. 4) are indicated by encirclement, (b) Equi- 
valent diagram notation for term -a 3 . 
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Figure F.9 - (a) Schematic notation whereby the elements of the determinant 
which contibute to equation (E.4) are indicated by encirclement, (b) Equi- 
valent diagram notation for term ^a . 
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Figure F. 10-(a) Schematic notation whereby the elements of the determinant 
which contibute to equation (E.4) are indicated by encirclement, (b) Equi- 
valent diagram notation for term ~a 3 
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Figure H. 1 - Diagram representation of I v . (a) Initially, 
(b) After the Fourier transform. 
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